Thermoelectric internal current loops inside inhomogeneous systems 



(N 

O 
(N 

a 

in 



Y. Apertetj^'B H. Ouerdane,^ C. Goupil,'' and Ph. Lecoeur'* 

^ Institut d'Electronique Fondamentale, Universite Paris-Sud, CNRS, UMR 8622, F-91405 Orsay, France 
'^CNRT Materiaux UMS CNRS 3318, 6 Boulevard Marechal Jmn, 14050 Caen Cedex, France 
^Laboratoire CRISMAT, UMR 6508 CNRS, ENSICAEN et Universite de Caen Basse Normandie, 
6 Boulevard Marechal Juin, 14050 Caen, France 
^Institut d'Electronique Fondamentale, Universite Paris Sud CNRS, 
91405 Orsay, France, CNRS, UMR 8622, 91405 Orsay, France 
(Dated: January 6, 2012) 

Considering a system composed of two different thermoelectric modules electrically and thermally 
connected in parallel, we demonstrate that the inhomogeneities of the thermoelectric properties of 
the materials may cause the appearance of an electrical current, which develops inside the system. 
We show that this current increases the effective thermal conductance of the whole system. We 
also discuss the significance of a recent finding concerning a reported new electrothermal effect in 
inhomogeneous bipolar semiconductors, in light of our results. 

PACS numbers; 72.20.Pa, 44.90.+C, 72.15.Jf 
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Thermoelectric power generation is a promising way 
to achieve efficient waste energy harvesting. To ensure 
a high heat-to-electrical power conversion efficiency, the 
thermal conductances of the materials used for thermo- 
electric modules (TEM) have, in principle, to be as low 
as possible[l|. Fu et al [1] recently reported on an elec- 
trothermal process that can modify the effective ther- 
mopower of semiconductor devices. In particular, they 
claimed that the joint application of a temperature gra- 
dient and an electric field (perpendicular to each other) 
to a bipolar semiconductor structure induces steady cur- 
rent vortices (even in open circuit configuration) which in 
turn yield Joule heating whose effect is to lower thermal 
conductivity. It is thus worthwhile to check whether this 
effect may be used to improve the so-called figure of merit 
ZT of bipolar semiconductor structures to a significant 
degree. 

The theoretical prediction of Fu et al i2] that inter- 
nal current vortices formed at a pn junction provide a 
way to reduce thermal conductivity in practical devices 
deserves closer inspection. In this Brief Report, using a 
macroscopic description of a two-leg TEM, we demon- 
strate that an internal current also gives rise to advective 
thermal transport which, unfortunately for practical ap- 
plications, largely compensates the effect proposed by Fu 
et al ^ and hence effectively lowers ZT. Studying the 
simple case of two thermoelectric modules connected in 
parallel both thermally and electrically, we suggest that 
internal currents caused by a temperature gradient are 
not directly linked to the transverse electrical field as 
supposed in Ref. but rather caused by thermoelectric 
inhomogeneities inside the materials. 

To gain insight into the main features of internal cur- 
rent loops, let us consider two thermoelectric modules 
TEMi and TEM2, and the equivalent module, denoted 
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FIG. 1. Coupled thermoelectric modules (top) and equivalent 
module (bottom). 



TEMoq, resulting from their association in parallel, both 
electrically and thermally as shown in Fig. [T] Each of 
them is characterized by its isothermal electrical conduc- 
tance Gi, its thermal conductance under open electri- 
cal circuit condition A'o,i and its Seebeck coefficient a;, 
where i can be 1,2 or eq as appropriate. All these co- 
efficients are supposed constant. The whole system is 
subjected to a temperature difference AT — Thot — ?coid, 
and its average temperature is T. 

Using linear response theory we can express each elec- 
trical current li and thermal flux Iq as functions of gen- 
eralized forces related to the temperature difference AT 
and voltage difference AV^ to which the TEM is sub- 
jected. The potential differences AV and AT are the 
same for TEMi and TEM2 and hence for TEMoq by con- 
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struction. The relation between the fluxes and the forces 
is given by Q : 
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from which we obtain a quite simple expression of the 
thermal flux: 



(2) 



This equation shows that two distinct processes con- 
tribute to the thermal transport: one is linked to thermal 
conduction by both phonons and electrons when there is 
no current flowing inside the structure (term in KqAT), 
the other to electrical current flow (term in aTI). Since 
this second contribution is associated with a macroscopic 
displacement of electrons it can be stated to be thermal 
transport by electronic advection; and the heat quantity 
transported by each electron ^] is given by |a|Te, e being 
the elementary electric charge. This notion of electronic 
advection is central to explain the increase of thermal 
conductance when an internal current develops inside the 
structure. We stress that the additional term should not 
be confused with the electronic part of Kq, which is used, 
for example, in the Wiedemann-Franz law. On Fig. [TJ 
this is shown with the added thermal conductance pa- 
rameterized by the electrical current. 

Besides constitutive laws for each module given by 
Eq. ([T]), there are additionnal relations linked to the par- 
allel configuration that must be accounted for. First, we 
ensure electrical current conservation; 
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Next, we consider that the mean thermal flux flowing 
through TEMeq is the sum of the two mean thermal 
fluxes flowing through TEMi and TEM2; hence 
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(4) 



Now using Eqs. (P), ^ and dU we are free choose the 
thermal and electrical configurations which permits an 
easy derivation of the equivalent parameters of the sys- 
tem as a whole. Under isothermal condition, AT — 0, 
and Eq. ^ leads to: 



Gi + G2, 



(5) 



Under closed circuit condition, AV = 0, and Eq. ^ 
leads to GcqCtcq = Giai + G2a2, so that the equivalent 
thermopower acq is defined as the weighted average of 
the two Seebeck coefficients ai and 



Giai + G2a2 
Gi + G2 



(6) 



This equation is the same as the one given by Hicks 
and Dresselhaus (Eq. (2) of Ref. for a semiconduc- 
tor with two conduction bands. The correspondence be- 
tween both is explained by the fact that each conduction 
band can be associated with a medium where electrons 
fiow parallelly. From this result we see that the effective 
Seebeck effect cannot exceed the larger one of the two 
materials. 

The equivalent thermal conductance is now determined 
under open circuit condition, /gq = 0. For nonzero values 
of Ii and I2, this condition is satisfied for Ii — ~l2', since 
the conservation of the thermal flux Q remains valid, we 
obtain: 
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(7) 



using Eq. ([2|). To proceed, we determine the intensity Ii 
as follows. Under open ciruit condition I^q = so that 
the voltage reads AV ~ — OfoqAT, which we include in 
the expression of Ii given by ([T]) to obtain: 
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This equation shows that there is a non-zero electrical 
current flowing as long as the two Seebeck coefficients 
are different. Now, the substitution of the obtained ex- 
pression of Ii into Eq. ^ , yields the equivalent thermal 
conductance at zero current: 
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(ai - a2fT (9) 



The above formula exibits an additional term next to 
the sum of the thermal conductance of each module. 
This term is related to the internal current that devel- 
ops inside the structure when sumitted to a temperature 
gradient. This current is proportional to the difference 
between the Seebeck coefficients of each leg. The total 
transported heat by this current is proportionnal to this 
difference too, so the increase in thermal conductance is 
proportionnal to (ai — 012)^. 

As an internal current is generated, energy dissipation 
is caused by Joule effect. The total dissipated power is 
given by the sum of the power dissipated in each part 
of the whole system: Pjouic — /^t(l/Gi -I- I/G2), where 
/int. (= ^1) is given by Eq. an explicit expression is. 
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(10) 



The dependence on the square of the difference of the 
Seebeck coefficients of the modules shows that the more 
they are dissimilar in terms of thermopower the more 
Joule dissipation is important. Here, there is no need to 
assume, as one would for the classical generator with two 
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legs, that each module has a different doping type. We 
also note that the internal current is stronger for mate- 
rials with high electrical conductance, hence metals are 
more sensitive to small inhomogeneities in the Seebeck 
coefficient. This explains why this effect is exploited in 
non-destructive testing to probe metallic inclusions in a 
host metal 

Let us remark that for inclusions at the nanoscale the 
previous statement about the increase of the thermal con- 
ductance no longer stands since these inclusions have a 
strong impact on thermal conduction by phonons. For 
example Kim et al demonstrated that nanoinclusions 
are efficient to lower the thermal conductivity in InGaAs. 

We now turn to the analysis of the paper of Fu et al [2] 
using our Eq. ([9|), which allows to explain the behavior 
of a pn junction submitted to a transverse thermal gra- 
dient under open circuit condition. A typical Pisarenko 
plot i.e. the Seebeck coefficient plotted against the 
carrier concentration, shows that the Seebeck coefficient 
is higher for low carrier density. In a pn junction a de- 
pleted zone forms and develops on each side of the in- 
terface over a few hundreds of nanometers depending on 
the doping concentration: in these two regions, the See- 
beck coefficients increase locally and thus become greater 
than those in the quasi-neutral regions. Since the carrier 
concentration inhomogeneity is transverse to the applied 
temperature gradient, this situation is similar to the sim- 
ple one studied in the present paper and we can expect 
internal current vortices to develop as described in 
We believe that the electric field present in the space 
charge zone is a consequence of the carrier depletion and 
is not a cause per se of the internal current generation: 
in a structure where one can tune the Seebeck coefficient 
without changing the carrier density we expect to find 
the same behavior. However, for the particular case of 
a pn junction the Seebeck inhomogeneity and the trans- 
verse electric field are closely connected through the car- 
rier depletion zone; therefore linking, in this case, the 
internal current either to thermoelectric inhomogeneties 
or to a transverse electric field essentially amounts to ex- 
press two viewpoints on the same phenomenon. Besides, 
the appearence of two vortices, one on each side of the 
junction, is due to the potentiel barrier arising at the in- 
terfaces: each type of carrier is confined in its own side 
so that these two separate systems have no infiuence on 
each other except when one of the sides becomes thinner 
than the depleted zone. 

We state further that a decrease of the thermal con- 
ductance is not possible in such circumstance: Fu and 
co-workers overlooked the advective part of the ther- 
mal fiux. According to these authors, the reduction is 
due to the fact that half of the energy dissipated by the 
Joule effect is actually going back to the hot side. This 
heat quantity sould be compared to the one transported 
by electronic heat advection: the ratio between these two 
quantities scales as AT/2T. In the framework of linear 
theory, the temperature difference AT is assumed to be 



small, so the heat flowing back to the hot reservoir is 
small compared to the advection part: internal currents 
increase the thermal conductance; they do not decrease 
it. Generation of internal currents in spatially inhomo- 
geneous systems such as multilayer structures were dis- 
cussed twenty years ago by Saleh et al [ll|; our analysis 
is consistent with their conclusions. 

The specific case of TEMs with both electrical and 
thermal parallel configurations satisfies Bergman's theo- 
rem, which states that the figure of merit ZT of a com- 
posite material cannot be greater than the larger ZT of 
the consituents [l^ • We have shown that the Seebeck co- 
efficient is lower than the larger one of the two, Eq. ([6|). 
The effective electrical conductance of two conductances 
in parallel is always greater than the larger one, but, 
since ZT cx G/Kq for a give temperature T, this in- 
crease is counterbalanced by the fact the thermal conduc- 
tances behave similarly when no internal current devel- 
ops. Moreover, considering a TEM designed with compo- 
nents that have different Seebeck coefficients, the increase 
in equivalent thermal conductance is even stronger. So, 
as expressed by Bergman and Levy p^ . the figure of 
merit of the TEMoq can only be lower than the highest 
ZT of the more efficient TEM. Indeed, the addition of 
an internal current can only lower the figure of merit, 
since it adds a dissipative process to the system. A full 
generalization of this result to composite materials ne- 
cessitates the derivation of an expression of the effective 
Seebeck coefficient for a configuration where TEMs are 
both thermally and electrically in series. Since a com- 
posite system may be viewed as a network of TEMs a 
general formulation can, in principle, be obtained; how- 
ever we anticipate that such derivation, which is beyond 
the scope of the present Brief Report, can be quite tricky. 
The simple example given in the present paper suffices 
to gain insight into the problem of internal currents and 
their properties. 

To end this Brief Report, we make two additional re- 
marks: 

Assuming that the properties of the thermoelectric 
modules are temperature-dependent, it would be inter- 
esting to check if it is possible to create instabilities in- 
side the module analogous to the Rayleigh-Benard con- 
vection phenomenon which occurs in conventional fiuids 
subjected to a thermal gradient. 

The present work and the cited ones on current 
loops bring us back to the original interpretation of the 
thermoelectric effect made by Seebeck who thought that 
the temperature difference led to magnetism whereas it 
was the internal current of the structure that created the 



magnetic field, which deviated the compass needle [13|. 
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